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Abstract 

The Jaynes-Cummings model, describing the interaction between a single two-level system and a 
photonic mode, has been used to describe a large variety of systems, ranging from cavity quantum 
electrodynamics, trapped ions, to superconducting qubits coupled to resonators. Recently there 
has been renewed interest in studying the quantum strong-coupling (QSC) regime, where states 
with photon number greater than one are excited. This regime has been recently achieved in 
semiconductor nanostructures, where a quantum dot is trapped in a planar microcavity. Here we 
study the quantum strong-coupling regime by calculating its photoluminescence (PL) properties 
under a pulsed excitation. We discuss the changes in the PL as the QSC regime is reached, which 
transitions between a peak around the cavity resonance to a doublet. We particularly examine the 
variations of the PL in the time domain, under regimes of short and long pulse times relative to 
the microcavity decay time. 
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Strong coupling is a phenomenon that has gained fundamental importance in the field 
of quantum technologies, where combinations of different systems are coupled quantum 
mechanically. Typically of such engineered systems, a matter field consisting of two levels is 
coupled to a photonic or phononic bosonic mode. For low densities, defined as being less than 
one excitation per two-level system, the new excitations of the system are superpositions 
of the matter and bosonic fields. The ability to create a superposition between matter and 
bosonic fields is important in terms of quantum information processing applications due to 
their ability to couple isolated quantum systems together. For example, a common bosonic 
mode may be used as a quantum communication channel between qubits, providing effective 
two-qubit interactions. The Jaynes-Cummings Hamiltonian, which describes any two-level 
system interacting with a single bosonic mode, is given by 

^sys = ^0 + ^int) (l) 

with 

■Ho = hua+a' + {huj- A) a^a, (2) 

"Hint = g {cr'a^ + cr^a) , (3) 

where cr^ are the raising and lowering operators for the two-level system, and is the 
creation operator for the bosonic mode. The coupling between the photon and the two level 
system is given by g, the transition frequency of the two-level system is u, and A is the 
detuning between the bosonic mode and the qubit jl-3|. 

At high densities, where there is more than one excitation per two-level system, the 
model predicts the quantum strong-coupling (QSC) regime. The quantum strong-coupling 
regime has been studied using superconducting qubits [4-l6|. This may be seen by directly 
diagonalizing the Hamiltonian, which may be easily performed since it preserves the total 
excitation number ritot = c^^a + |-^)(-^|- Considering the case of a resonant (A = 0) qubit 
and cavity system, the eigenstates may be found to be 

\±,n) = ^{\X)\n-l)±\G)\n)) (4) 

where \n) = '^{(^^)"'\^) ■ The ground and excited states of the two-level system are denoted 
\G) and \X), respectively. The energies of the states are 

= nhu ± 2g^/n (5) 
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thus for each excitation number n (with the exception of n = which is unique) there are two 
solutions forming the Jaynes-Cummings ladder, with the energy splitting increasing with the 
excitation number n. Recently an experiment studying the Jaynes-Cummings ladder in the 
quantum strong-couplingregime was realized, where a single quantum dot was placed in a 



microcavity structure 



8| . The highly excited states of the Jaynes-Cummings ladder were 
accessed by a combination of off-resonant pumping which first excites the ntot = 1 manifold, 
then resonant pumping which results in climbing to higher ntot states. 

In this paper we perform a numerical study of the photoluminescence of the quantum 
dot-microcavity system. While similar studies of the photoluminescence have been made 
for the Jaynes-Cummings model in a number of works js ili, the pumping scheme that is 
used in these works typically consider a continuous-wave pumping. However, in order to 
achieve the high densities required for the QSC regime, experiments typically use a pulsed 
excitation due to the power constraints of the typical lasers used. Such pulsed excitations 
give a time- dependent, non-equilibrium element to the problem. The dynamics of the system 
is determined not only by the Hamiltonian ([1]), but due to the open-dissipative nature of 
the matter and bosonic fields, which have a finite decay time. In our calculations, we 
consider both regimes where the pulse duration is shorter and longer than the lifetime of the 
particles in the system. Due to the pulsed excitation, the system is in general always in a 
non-equilibrium state. This allows for the possibility of directly studying the system during 
the decay process via observation of the PL characteristics. 

We also consider the effects of various pumping configurations, namely with respect to 



either the matter or bosonic fields 
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13| . In general, it is also possible to observe the 



photoluminescence that is associated with either the matter or bosonic fields. In the case 
of experiments in Ref. [3], either the PL due to the photons leaking out of the cavity, or 
the PL due to recombination of the excitons can be distinguished. These can in principle 
give different spectral characteristics, as we find explicitly from our calculations. This gives 
four possible pumping/measurement schemes, the relevant one depends on the particular 
experimental configuration used. We calculate all four cases and find that there are particular 
spectral characteristics common to all configurations which give spectral evidence for the 
QSC regime being reached. 
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FIG. 1: The particles, oscillating between the ground state \G) and the excited state \X), are 
pumped in the cavity by a intense pump laser, either Pq or Px- The excited particles in the cavity 
decay with either the exciton decay rate 7 or cavity photon decay rate n. The coupling strength 
between the cavity photon and a two-level atom is given by g. 



RESULTS 



To model the open-dissipative dynamics of the coupled cavity-quantum dot system, we 
use a master equation with the coherent part obeying ([T]), Lindblad photon and exciton 
decay terms with rates k, and 7 respectively, and Lindblad photon and exciton pump terms 
with rates Pc and Px respectively (see Figure [1]). The pump terms are assumed to obey 
a Gaussian profile in the time domain with time duration Tp (see Methods). The density 
matrix elements are numerically evolved in time, and the photoluminescence is calculated 
using two-time correlation functions evaluated using the quantum regression theorem jo- 
A cutoff rzinax foT the number of photons is made sufficiently large such that the trace of the 
density matrix is close to unity for all times. For concreteness we have chosen parameters 
that are of the order of the experiment performed in Ref. Tj. Specifically, the coupling 
constant g is of the order of 10 meV, while the excitation energy hw is of the order of eV; 
thus we fix (7 = 10 and hw = 1000 where all the parameters are in units of meV. The 
detuning parameter is set to A = throughout this study. 
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FIG. 2: Density matrix dynamics p{t) for long pulse pumping with Tp > 1/k. The time evolution of: 
(a) the occupation PGn,Gn for the ground state Gn under photon pumping Pc, (b) the occupation 
PXn,Xn for the excited state Xn under photon pumping Pq, (c) the occupation PGn,Gn for the 
ground state Gn under exciton pumping Px, and (d) the occupation pxn,Xn for the excited state 
Xn under exciton pumping Px- Curves are plotted in order of red (n = 0), green (n = 1), blue 
(n = 2), red (n = 3), green (n = 4), blue (n = 5), etc., with photon number n. As references, the 
dashed curves indicate the pumping profile. Parameters used here are g = IQ, Ttw = 1000, Tp = 10, 
7 = 0, and k = 1. Negative times refer to prior to the pulse at t = 0. 

Density matrix dynamics 

We consider two regimes in this study, where the pump timescale Tp is either much longer 
or shorter than the polariton lifetime I/k, when the system is pumped by either Pc or Px- 
Considering the case of long pumping (rp > 1/ k) first, we show the time evolution of the 
populations and the pump-pulse profile in Figure [2J The pumping pulse is centered at to = 0, 
and the system starts in the vacuum state \G) |n = 0). We see that as the pump is increased, 
the probability that the state is in the vacuum state decreases, resulting in excitations being 
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FIG. 3: Density matrix dynamics p{t) for short pumping pulse with Tp < 1/k. The time evolution 
of: (a) the occupation PGn,Gn for the ground state Gn under photon pumping Pq, (b) the occupa- 
tion pxn,Xn for the excited state Xn under photon pumping Pq, (c) the occupation PGn,Gn for the 
ground state Gn under exciton pumping Px, and (d) the occupation pxn.Xn for the excited state 
Xn under exciton pumping Px- Curves are plotted in order of red (n = 0), green (n = 1), blue 
(n = 2), red (n = 3), green (n = 4), blue (n = 5), etc., with photon number n. As references, the 
dashed curves indicate the pumping profile. Unlike the case with long pulse, the system evolves 
asymmetrically with time. Parameters used here are g = 10, hco = 1000, Tp = 1.5, 7 = 0, and 
K = 0.1. Negative times refer to prior to the pulse at t = 0. 

created. Cutoffs of n^ax = 50 for Pg and n^aax = 30 for Px are sufficient for the density 
matrix to have Tip = 1 for all times. At the strongest pumping intensity, there is a dip in 
some of the probabilities (e.g. the |X)|n = 1) state), signifying that excitations beyond the 
ntot = 1 manifold are being excited. This is by definition the QSC regime, and confirms 
that with suitably intense pumping such regimes may be reached. 
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Similarly, Figure |3] shows the density matrix elements for the case of a short pulse time 
Tp < 1/k,. Again we see that many of the states above the ntot = 1 manifold are being 
excited, confirming that either pumping scheme may be used to reach the QSC regime. We 
however see an asymmetric behavior to the density matrix elements, unlike Figure [2] which 
was approximately symmetric. This can be explained due to the relative timescales of the 
pumping, where either the system can respond fast enough to changes in the pumping (the 
former case), and the pump is faster than the timescale of the system (the latter case). In 
the case of the short pulse, the polaritons require a time ~ 1/k before they can escape the 
system, thus the decay process itself is being observed for times t > to. 

Photoluminescence spectra 

We now show the photoluminescence dynamics of the quantum dot system. We consider 
two cases of long pump pulses (rp > I/7) and short pump pulses (rp < 1/k, I/7). For 
each case we evaluate the four cases of cavity photon pumping and exciton pumping with 
the measurement of the photonic or excitonic PL (Figure [Hand [6]). From the evolved matrix 
elements we also calculate the corresponding expectation values for the photon number 
^ph = o^a and the exciton number = (Figure [5] and [7]). We assume that the 

exciton decay rate 7 is negligible when measuring the photon mode, whereas the cavity 
photon decay rate k is negligible when measuring the exciton mode, such that 7 = and 
K = for photonic PL and excitonic PL, respectively. 

Photon PL and photon pumping 

First we examine the case where the system is photon pumped, and decays via photon 
emission k = 1. The average photon and exciton number is shown in Figure [5]( a) for the 
long pump pulse case. We see that the photon number is much greater than the exciton 
number at the peak densities, which is a consequence of there being a maximum exciton 
number of 1, whereas there is no limit to the maximum photon number. For the quantum 
dot system this is a statement of Pauli's exclusion principle, as two excitons may not occupy 
the same dot (assuming the dot size is the same or smaller than the exciton Bohr radius). 
At high densities, there are therefore many more photons than excitons, again illustrating 
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FIG. 4: Photoluminescence spectra for long pump pulses. The time evolution of: (a) photonic 
PL by photon pumping Pc, (b) photonic PL by exciton pumping Px, (c) excitonic PL by photon 
pumping Pc, and (d) excitonic PL by exciton pumping Px in the case of long pump pulse. The 
peak PL approaches the cavity photon mode at high densities, whereas two peaks at UP and LP 
are brighter at low densities. Dashed line indicates the cavity photon energy. Parameters used are 
g = 10, huj = 1000, TP = 10, and: (a) P^'^'^ = 0.8, «; = 1, 7 = 0, n^ax = 50; (b) P^^^ = 6.0, k = 1, 
7 = 0, nn,ax = 30; (c) P^^^ = 0.3, k = 0, 7 = 1, n^a.. = 120; and (d) P^^^ = 3.5, k = 0, 7 = 1, 

'^max — 40. 

the QSC regime. For very high densities, such that (riph) ^ 1, the bosonic operators of the 
Hamiltonian ([1]) may be written as a c-number, giving the effective Hamiltonian [l^ 



huj{ntot) +9\/{nph) cr' 



(6) 



which has eigenstates ^ ± \G)). Thus at high density the exciton population should 
be at most (riex) = 1/2. In Figured we see that the exciton number never exceeds this 
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FIG. 5: The mean number of cavity photons (up^) and excitons (ncx) for long pump pulses. The 
cases shown are (a) photonic PL by photon pumping Pc, (b) photonic PL by exciton pumping 
Px, (c) excitonic PL by photon pumping Pc, and (d) excitonic PL by exciton pumping Px in the 
case of long pump pulse. Parameters used are g = 10, hco = 1000, rp = 10, and: (a) p™^^ = 0.8, 
K = 1, ^ = 0, nmax = 50; (b) P^^" = 6.0, k = 1, 7 = 0, n^nax = 30; (c) P^^" = 0.3, k = 0, j = 1, 
n^nax = 120; and (d) Pj^^" = 3.5, k = 0, 7 = 1, n^ax = 40. 

value, in agreement with this argument. 

In Figure Hl^a) we see that at the highest densities the brightest peak comes from the cavity 
photon energy. As the system decays via leakage through the microcavity, the peak moves 
towards the lower polariton (LP) and upper polariton (UP) resonances, finally resulting in 
only the LP and UP resonances at the lowest density. A qualitative understanding of this 
spectrum can be described by using the dressed states |±,n) defined in equation (j4]). The 
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strong photon pumping initially creates a state which is 

(at)"|G)|0) oc 1^)1^) 

= ^(|+,r^)-|-,r^)). (7) 

As the photon PL is observed due to a loss of a photon out of the system, this state transitions 
to 

a\+,n) ^ y/n |+,n - 1), 

a\ — ,n)^^/n\ — ,n — l), (8) 

where we have assumed n ^ 1. The transition energy is then 

AE = - Eti ^ nw, (9) 

which corresponds to the cavity photon energy. In the low density limit, we have energy 
peaks at UP and LP as shown in Figure HJ^a). Thus, we find that the energy at the cavity 
photon arises from the enhanced transitions of 1+,"^) — )■ \+,n — 1) and \ — ,n) — t- \—,n — 1) 
due to the increased number of photons. The initial state has the same amplitude for the 
states \+,n) and \—,n), which is the reason for the symmetrical energy peaks to the cavity 
photon mode. 

In Figure [7] we plot the average number of photons and excitons for the short pump pulse 
case. We again see an asymmetric distribution of the particle numbers with respect to the 
maximum pump intensity, due to the slower decay time of the system with respect to the 
pump time. In all cases we again observe that the exciton number never exceeds 0.5 due to 
the same reasons as the long pump case. 

In Figure |6] we show the PL spectrum for short pump pulses. We again see that at the 
highest densities the dominant peak is around the cavity photon energy. As the system 
decays away, there is small shift of the peak in the direction of the LP or UP energy. At 
the lowest densities, only the LP or UP peak is visible. Apart from the decay dynamics, the 
PL signatures are fairly similar for both the long and short pump cases. In both cases we 
note that despite the high average densities achieved as seen in Figure Ej there is always a 
resonance peak visible at the LP and UP energies. According to the simplistic treatment in 
equation ([8]) we would expect that these peaks would only be seen at the lowest densities. 
We see that in fact the dynamics is more complicated than the simplistic picture by again 
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considering Figures |2] and [31 If it were purely the case of ([8]) we would only see the density 
matrix elements corresponding to the highly excited states and zero population of the low 
excitation states such as Pgi,gi and pxo,xo- In fact we see that for both the long and short 
pump-pulse cases, there is no reordering of the states, so that low excitation states are always 
more probable than the higher-excitation states. Thus there is always a strong contribution 
of the transitions |±, 1) |0) which are simply the LP and UP transitions. For this reason, 
for the typical pumping strengths that we consider, the PL spectrum consists of a mixture 
of both high and low excited states. Only the parts of the PL spectrum in the vicinity of 
the cavity photon energy correspond to states in the QSC regime. 

Photon PL and exciton pumping 

We now discuss the case of exciton pumping Px with photon decay, where we fix the 
decay rates to k = 1 and 7 = 0. The time evolution of the PL is shown in Figure Hl^b) 
and the populations in Figure O^b). In comparison to the photon pumping case, the system 
requires high pump powers to reach the high density due to the spin operator suffering 
from phase-space filling, {a'^)"'\0) = {n > 2). In order for the system to reach high densities 
for this pumping configuration, the pump must first create an exciton in the system, then 
this must quickly become transfered to a photon via the coupling between the cavity and the 
quantum dot. If the pump can create another exciton in the system before the first photon 
decays, then high densities may be reached. The condition required for reaching the QSC 
regime is therefore Px,g ^ i^- Apart from this difference, we see a qualitatively similar 
spectrum to the photon PL/photon pumping case at high density, with the PL emerging 
most strongly at the cavity photon resonance. As the system decays, there is a shift of this 
peak towards the LP and UP resonances. 

For the short pump pulse case, we see an almost identical spectrum to the long pump-pulse 
case [Figures int^b) and[7l^b)]. We can therefore conclude that despite the more inefficient 
pumping method of the exciton pumping, as long as the criterion Px, g ^ k is satisfied, the 
PL and population characteristics are independent of the particular pumping scheme that 
is used. 
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Exciton PL and photon pumping 



We now turn to the excitonic PL and photon pumping case and fix the decay rate k = 
and 7 = 1. As seen in Figure life), we observe a quahtatively different spectrum, with 
stronger side peaks at the LP and UP resonances, as well as the cavity photon resonance. 
We understand this by considering our simplified picture again. For strong photon pumping, 
we start once more in the state defined in equation ([7]), {a'^)"'\G)\0) oc (1+,'^') — | — 5^^))- 
As the exciton PL corresponds to the loss of an exciton, we consider the transitions induced 
by operating a~ to the above states 



Unlike the photon PL case in equation ([8]) where \^,n) only transitions to one state, here 
there are two possible states to jump to. In this case all the transitions between \^,n) — )■ 
|±,ra — 1) have the same amplitude. Therefore, the most distinctive feature here is that 
we have side peaks at AE ^ hu ± 2g^fn^ in addition to IS.E ^ hw for high density. The 
kind of transitions shown above are familiar from resonance fiuorescence, where in the high- 
excitation case the photoluminescence spectrum reduces to the Mollow's triplet spectrum 

[xsj. The difference here is that the densities are not quite so high, such that we may treat 
the excitation field as a classical field, and also the non-equilibrium nature of the pumping 
profile. Nonetheless, similar spectral characteristics may be seen for the high-excitation 
regime we consider here. For continuous and high-pumping rates we have confirmed that 
our numerics reproduces the Mollow's triplet spectrum. 

For short pump pulses, as shown in Figures [6](c) and O^c), the numerical simulations 
are less stable with respect to the pumping dynamics, hence we were only able to achieve 
relatively low densities (although still in the QSC regime) in comparison to the long pump- 
pulse cases. Nevertheless, we still see the characteristic Mollow's triplet spectrum evolving 
towards the LP and UP resonances. From this we conclude that both the short and long 
pumping profiles should be effective in reaching the QSC regime, giving similar spectral 
characteristics. 




(10) 
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Exciton PL and exciton pumping 

We finally consider the case where excitons are pumped and the exciton PL is measured 
with 7 = 1 and k = 0. The time dependence of the excitonic PL is shown in Figures HJ^d) 
andE](d). We see a qualitatively similar PL spectrum to the excitonic PL and cavity photon 
pumping case, with prominent side peaks on either side of the cavity photon resonance. 
As was the case with photonic PL with excitonic pumping, the system requires high pump 
powers in order to reach the high-density regime, and requires Px,g ^7- In a similar way 
to the exciton PL and photon pumping case, all the four transitions in equation ffTOj) occur, 
giving rise to more prominent side peaks. We thus conclude that the particular form of the 
pumping from either short or long pulses, photonic or excitonic pumping gives qualitatively 
similar PL spectra for the excitonic PL cases. 

DISCUSSION 

We have obtained the PL spectrum of a single quantum dot in a cavity pumped via 
injection of photons or excitons in a pulsed Gaussian profile. The particles decay either 
with the cavity photon decay rate k, or with the exciton decay rate 7 depending on how 
we measure the PL, i.e. photonic PL or excitonic PL. We first find that with any choice of 
pumping scheme the quantum strong-coupling regime can be reached for sufficiently strong 
pumping strengths. While this is intuitively clear for the photonic pumping case, this is 
also true for the exciton pumping, which suffers from the Pauli exclusion principle excluding 
it from the simultaneous injection of more than one exciton at a time. Despite this, the 
high density regime may be reached as long as the time required for conversion of excitons 
to photons via the coupling is faster than the decay rate P,g^K,, 7. We also generally 
find that the PL spectrum that is measured is fairly insensitive to the particular pumping 
method used to excite the system. Our considerations of exciton and photon pumping, 
using both short and long profiles gives fairly consistent spectra throughout. Naturally, if 
the decay dynamics itself needs to be probed, then using a short pulse is really the only 
option. However, if the only goal is to reach the QSC regime, then either case can be used. 

The general trend of the photonic PL is that there is a peak at the cavity photon mode at 
high densities and as the system decays away it moves towards the UP and LP resonances 
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FIG. 6: Photoluminescence spectra for short pump pulses. The time evolution of: (a) photonic 
PL by photon pumping Pc, (b) photonic PL by exciton pumping Px, (c) excitonic PL by photon 
pumping Pc, and (d) excitonic PL by exciton pumping Px- The peak PL approaches to the cavity 
photon mode at high densities, whereas two peaks at UP and LP are bright in addition to the 
peaks near the cavity photon resonance at low densities. Dashed line indicates the cavity photon 
energy. Parameters used are g = 10, huj = 1000, rp = 1.5, and: (a) p™^^ = 0.45, k = 0.1, 7 = 0, 
nmax = 40; (b) Pjf^^ = 1.3, K = 0.1, 7 = 0, n^,, = 30; (c) P^^^ = 0.45, k = 0, 7 = 0.1, n^^x = 40; 
and (d) Pf^'' = 1.3, K = 0, 7 = 0.1, n^^ = 40. 

in the low-density regime. For excitonic PL, the Mollow-triplet-hke spectrum occurs as 
the density is increased. These properties agree with the theoretical predictions of [l^ at 
high density. Some differences to the PL spectrum is seen at intermediate densities (Figure 
3(b) of Ref. [1^). In contrast to the many-exciton model, where the photoluminescence 
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FIG. 7: The mean number of cavity photons {rip^) and excitons (riox) for short pump pulses. The 
cases shown are (a) photonic PL and photon pumping Pq, (b) photonic PL and exciton pumping 
Px, (c) excitonic PL and photon pumping Pq, and (d) excitonic PL and exciton pumping Px- 
Parameters used are g = 10, fkj = 1000, rp = 1.5, and: (a) P™^^ = 0.45, k = 0.1, 7 = 0, 
nmax = 40; (b) P^^^ = 1.3, k = 0.1, 7 = 0, n^ax = 30; (c) P^^^ = 0.45, k = 0, 7 = 0.1, n^^x = 40; 
and (d) P^^^ = 1.3, K = 0, 7 = 0.1, n^^ax = 40. 

smoothly transitions from the LP to cavity photon energy, in the single quantum dot case 
the PL discontinuously shifts between the two limits. We have confirmed that by increasing 
the number of excitons in the system the transition gradually becomes smoother jl^]. The 
discontinuous spectrum can therefore be seen as a signature of the quantum dot system, 
although in a typical experimental situation it is more likely that many quantum dots are 
excited simultaneously in order to have sufficient signal amplitude in the measurement. 
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METHODS 



Time evolution of the density matrix 

Due to the open-dissipative nature of the quantum dot-microcavity system, the PL char- 
acteristics are modeled by a master equation following Refs. js - 11|. We consider the reduced 
density matrix 

= [p, "Hsys] + 2 (Sapa'f - a)ap - pa^a) 

+— {2a~ pa^ — a^a~ p — pa^a^) 
P (t) 

H — (2a^pa — aa) p — paa)^ 

+ (2a+pa- - a-a+p - pa'a^) , (11) 

where we assume that the pump has a Gaussian time profile 

Pc{t), Px{t) = Po exp (-^^^) • (12) 

Using the basis \in) where i = G is the state with no exciton and i = X is the state with 
an exciton, we may define the matrix elements of the reduced density matrix 

Pinjm = {in\p\jm), (13) 

where n, m is the photon number. The time evolution of the density matrix elements are 
d i ^ 

-^PGn,Gn — J^QV n [pGn,Xn-l — PXn-l,Gn) + 1 PXn,Xn 
-K [npGn,Gn -(« + !) PGn+l,Gn+l\ 
— PxPGn,Gn + PcnpGn~l,Gn~l 

-Pcin + l)pG„,Gn, (14) 

d i . 

-^PXn,Xn = '^dVn + 1 {pXn,Gn+l — PGn+l,Xn) — lPXn,Xn 
-K [npxn,Xn - {n + 1) PXn+l,Xn+l] 
+PxPGn,Gn + PcnpXn-l,Xn~l 

-Pc{n+ l)pxn,Xn, (15) 
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-^PGn,Xn-l — -^[gV^T' iPGn,Gn — PXn-l,Xn-l) 
+ ^PGn,Xn-l] 
- [(7 + H{2n - 1) + Px) /2] PGn,Xn-l 



+K^/n{n + l)pGn+l, 



Xn 



+Pc 

y n{n — \)pQri-l,Xn—2 

Pg 

—{2n + l)pGn,Xn-l, (16) 



-rPXn-lfin — —-z\gV^iPGn,Gn — PXn-l.Xn-l) 

at n ^ 



+'^PXn-l,Gn] 
- [(7 + K{2n - 1) + Px) /2] pXn-l,Gn 
+ K^Jn{n + l)pXn,Gn+l 

+Pc^n{n - l)p 



Xn-2,Gn-l 



Pc, 



-^{2n+l)pxn-l,Gn. (17) 

Imposing a maximum photon number rimax we may evolve this set of equations which form 
a set of closed ~ 4nmax equations. In order to investigate the system dynamics, we truncate 
up to finite times instead of calculating up to the stationary limit. 

Two-time correlation function 

In the quantum dot system, two types of photoluminescence mechanism can be consid- 
ered: the particles are supplied from the pump laser and they decay due to leakage through 
the microcavity mirrors, or via exciton recombination [l^. We call these processes the 
photonic and excitonic PL, respectively. The PL is given by [l| 

S{r,v) = -Re / dr e'"^ G^^\t,T). (18) 
Jo 

where G^^\t, t) is the two-time correlation function, which for the quantum dot system, is 
either 

Gg^(t,r) = {a\t + T)a{t)), (19) 
G«(t,r) = (a+(t + r)a-(t)), (20) 

17 



depending on which type of PL is being calculated. The Fourier transformation of G^(p{t, r) 



gives the photonic PL which is emitted by the leakage from the imperfect mirror, while that 
of G'i\t,T) gives t 
Following Refs. 
operators 



le excitonic PL which is emitted due to the recombination of excitons. 
9l-lll|. in the interaction picture we may conveniently define the following 



\Xn){Gn\e'^^/\ 

\Gn + l){Xn-l\e'^^-^^^^/''. 



(21) 



In terms of the new operators defined in equation (ETJ), we can rewrite the two-time corre- 
lation function as 



G^c\t,r) = ^v/^(^(at.Jt + r)a(t)) + (at,„(t + r)a(t)) 



G'i\t,r) = 5^(a+(t + r)a-(t)). 



(22) 
(23) 



In order to calculate the two-time correlation function, we need to obtain the expectation 
values of two operators at different times. The time evolution of the expectation values for 
each operators is then 



|:(4n(r)) = (a^Jr)) 



h 2 



+ (4n+lM)'^V(^+l)(^ + 2) 



+ 



t 



(«Ln-2(r))J Pc^n{n-l) 
(«L-iM)]'^(2^ + l), 



(24) 
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dr 



lUJ — 



+ (a+_i(r))Pcri-(a+(r))Pc(n + l), 



(25) 



(4n-iW)^x-(a,^(r))-^v^ 



+K„-iM) 



^(;k.-A)-7-^(2n-l) 



+ (a^„(r))fi:v/^(n + 1) + (?„(r))-(?Vn + 1 
+ (4n-2M)^cVri(n-l) 



-(«k-lM)^(2r^ + l), 



(26) 



— — A Kn 

fi ' 2 



+ (^-n+i (r ) ) K A/ra(ra + 2) 



+ {<;n^^{T))Pc^{n-l){n + l) 
-{<,n{r))Pc{n + l). 



(27) 



Applying the quantum regression formula to equations (Eljl - flTTj) . we can compute and obtain 
the following sets of two-time correlation functions 16]. For the calculation of G^^\t, r), we 
need to obtain 



{(af;„(t + r)a(t)),(a+(t + r)a(t)), 



(28) 



and, for the calculation of Gx\t, t 



{(aL(t + r)(7-(t)),(a+(t + r)a-(t)), 



.i(t + r)(T-(t)),(.;„(t + r)a-(t)) 



(29) 
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is needed. The above equations satisfy the initial conditions 



Vn+ lpGn+l,Gn+l {t) , 



{a 



'Xn-l 



(t)a(t)) 



V^PGn+l,Xn(i)e'^*/", 



(30) 



and 



PXn,Gn+l 





:n-i(tMt)) 

(^n(t)a(t)) 



= 0, 




(31) 



for the calculations of GQ{t,T) and Gy(t,T), respectively. Here, r is calculated up to the 
stationary limit so that we can compute the two-time correlation functions. 
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